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I.  INTRODUCTION 

In  this  report,  a  method  for  computing  the  electromagnetic  scattering 
properties  of  complex  objects  is  presented.  This  method  can  deal  with 
electrically  small  or  large  bodies  alike  and  is  not  restricted  to  perfect 
conductors.  It  is  an  iteration  method  based  on  an  integral  representation  of 
Maxwel 1 1 s  equations. 

In  Chapter  2,  a  brief  review  of  the  integral  equations  for  the 
electromagnetic  field  is  given.  It  is  shown  how  they  can  be  combined  with 
impedance  boundary  conditions  in  order  to  deal  with  non-perfect  conductors  and 
as  a  special  case  of  interest,  coated  perfect  conductors. 

Chapter  3  takes  a  look  at  the  AS-MM  hybrid  diffraction  technique  which  is 
the  starting  point  of  the  iteration  method.  In  this  hybrid  technique,  the 
current,  induced  on  the  surface  of  a  scatterer  by  an  incident  field,  is 
approximated  by  optics  currents  which  are  supplemented  by  correction  currents 
in  order  to  give  the  true  current.  The  approximate  optics  currents  are 
computed  by  iteration,  whereas,  the  correction  currents  are  computed  by  using 
the  moment  method  in  certain  "moment  method  regions".  Although  this  approach 
has  been  shown  to  be  successful  for  simple  geometries,  the  task  of  dividing  up 
a  complex  body  into  moment  method  and  asymptotic  regions  is  a  formidable 
one.  These  moment  method  regions  complicate  the  programming  even  for  perfect 
conductors  let  alone  non-perfect  ones. 

In  an  attempt  to  do  away  with  the  moment  method  regions,  the  iteration 
method  was  developed.  This  method  is  described  in  Chapter  4.  The  two  methods 
are  basically  similar  in  that  they  both  define  the  same  approximate  optics 
currents  and  the  same  correction  currents.  However,  no  moment  method  regions 
are  defined  and  the  correction  currents  are  found  by  iteration  just  as  are  the 
optics  currents. 


1 


.>v 


In  Chapter  5,  some  numerical  results  for  the  two-dimensional  problem  of  a 
square  cylinder  are  exhibited.  Results  for  the  iteration  method  are  compared 
to  results  obtained  by  the  method  of  moments,  which  is  considered  exact.  The 
results  are  in  general  excellent  except  for  small  regions  in  the  shadow  side 
for  large  bodies.  This  question  is  examined  in  Chapter  6  and  a  possible  way 
of  improving  the  iteration  method,  by  changing  the  iteration  procedure  and  the 
approximations  made,  is  suggested.  Chapter  7  is  devoted  to  a  summary  and  a 
few  concluding  remarks. 


2.  THE  INTEGRAL  EQUATIONS  FOR  SCATTERING  FROM  NON-PERFECT  CONDUCTORS 


Consider  an  arbitrarily  shaped  non-perfect ly  conducting  body  situated  in 
a  vacuum.  Let  the  body  be  illuminated  by  an  incident  electromagnetic  field 
(f1,  ft1  ).  Then  the  total  electric  and  magnetic  fields  at  an  observation 
point  ft  are  given  by  (with  suppressed  time  variation  exp( jwt) )[1] 


£(ft)  =  Tft^ft)  -  T  {  {J8n0(n‘xfi)  G(r)  -  (n'x^xV*  G(r) 

s 

-  (n1  •  t)  V 1  G(r)>  ds'  (1) 

fyft)  =  TfV(ft)  +  T  ^  {j  —  (n'xft)  G(r)  +  (n’xft)  xV‘  G(r) 
s  no 


+  (n1  •  ft)  V'  G(r)}  ds' 


(2) 


where  T  =  (1  -  ft/4ir)"*  and  i  is  used  to  denote  the  principle  value  integral 

over  s,  i.e.,  the  integral  over  the  closed  surface  excluding  an  e  neighborhood 
of  the  singularity  in  G(r)  at  r  =  0.  The  absolute  value  of  the  solid 

angle,  ft,  subtended  by  the  deformed  surface  (in  taking  f)  at  the  field 

point  ft  must  be  determined  from  the  geometry  of  the  surface.  However,  if  s  is 
smoooth  at  ft  then  ft  =  2v  (T  =  2)  and  if  ft  is  not  located  on  the  surface 
then  ft  =  0  (T  =  1).  All  the  parameters  and  variables  (8,  nQ,  etc.)  connected 
with  the  electromagnetic  field  are  defined  in  Appendix  A.  Note  that  ft  is 
allowed  to  be  located  on  a  surface  whose  tangent  is  not  a  differentiable 

function  of  position,  however,  the  field  is  required  to  possess  a  finite  mean 
value  at  ft.  Equations  (1)  and  (2)  are  referred  to  as  the  electric  field 

integral  equation  (EFIE)  and  the  magnetic  field  integral  equation  (MFIE). 
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Introducing  the  equivalent  sources: 


3  =  nxfi,  A  =  -nxi:  (3) 


Equations  (1)  and  (2)  can  be  rewritten  as 

£(ft)  =  T£7(ft)  -  T  {  U3n03  G(r)  +  ftx§'  G(r)  -  (n‘  •  t)  G( r ) }  ds*  (4) 

s 

fi(ft)  =  Tfi1 (ft)  +  T  f  —  ft  G(r)  +  3x§'  G(r)  +  (n‘  •  fi)  G(r)}  ds‘  (5) 
s  no 

where,  from  Maxwell's  equations: 

n  .  t  =  ^  •  3,  n*  •  ft  =  J-  $•  •  ft  (6) 

o  pnQ 

The  scattered  field  can  be  obtained  from  Equations  (4)  and  (5)  as  follows: 


£s(&) 

=  £(&) 

-  £’(R) 

(7) 

-  fi(ft) 

-  fi^R) 

(a) 

and  since 

for  $ 

t  s,  T  »  1, 

Equations  (4)  and 

(5) 

together  with  (7) 

and  (8) 

yield 

£s(6) 

_  f 

~  ~  J 

s 

{j3nQ3  G(r) 

+  ftx$' 

G(r)  -  (n1 

•  £) 

V'  G( r) }  ds' 

(9) 

fis(S) 

-  /  (- 
s 

ft  G(r )  + 
no 

3xv' 

G{r)  +  (n1  • 

fi)  ? 

G( r ) }  ds' 

(10) 
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The  equivalent  surface  currents  can  be  found  from  Equations  (1)  and  (2)  by 
considering  $  e  s  with  s  smooth,  so  that  T  =  2  and  a  unit  vector  n  can  be 
defined  at  $  on  s.  Then  upon  taking  the  vector  cross  product  of  n  with  (1) 
and  (2)  and  using  Equation  (3),  Equations  (1)  and  (2)  yield: 


-fi(fi)  =  2  nx!1^)  -  2  nx-j-  {j6n0^G(r)  +  ftxV'G(r)  -  (n‘  •  £)  V'G(r)}  ds'  (11) 

s 

3(8)  =  2  nxfi1  (R)  +  2  nx  j{-  &  fiG(r).  +  3xV'G(r)  +  (n*  •  ft)  V'G(r)}  ds'  (12) 

s  no 

with  n'  •  £  and  n*  •  ft  given  by  Equation  (6). 

Equations  (11)  and  (12)  are  two  coupled  equations  for  the  surface 
currents  3  and  ft,  which  when  computed  can  be  used  in  Equations  (9)  and  (10)  to 
find  the  scattered  field.  These  equations  can  be  decoupled  by  using  the 
impedance  boundary  conditions  [2]: 


-is 


On  substituting  (13)  and  (14)  in  (11),  (12),  and  (6),  the  following  decoupled 
equations  are  obtained: 


-fi(ft)  =  2nx£1  ($)  -  2nx  f  {^(nxft)  G(r)  +  ftxv*  G(r)  -  (n1  •  t)v'  G(r)}  ds'  (15) 

s  r 

3(8)  =  2nxft1 (ft)  +  2nx  ^{j0Zp(nx3)  G(r)  +  3x$'  G(r)  +  (n*  •  ft)  V*  G(r)}  ds' (16) 


ith  n'  •  £  and  n'  •  ft  given  by: 


(17) 


n1  •  £  =  $‘.(n'xft),  n-fi  =  $'-(n'x3) 

Where,  Zp  =  Z/nQ  is  the  relative  surface  impedance  of  the  scatterer. 

The  complexity  of  the  above  integral  equations  can  be  reduced  by 
considering  a  specific  physical  situation.  The  scatterer  that  will  be 
considered  here  is  that  of  a  perfect  conductor  coated  with  a  layer  of 
non-perfectly  conducting  material.  Furthermore,  the  following  assumptions 
will  be  made  [3] : 

1.  The  index  of  refraction  of  the  coating  is  very  large  and  has  a  large 
imaginary  part.  |N|  >>  1  and  Im(N)  >>  1. 

2.  The  fields  outside  the  body  vary  slowly  with  respect  to  the  dimension  of 
the  interior  wavelength. 

3.  The  thickness,  d,  of  the  coating,  is  less  than  the  principal  radii  of 
curvature  of  the  outer  surface,  and  8d  <<  |N| 

4.  The  radii  of  curvature  of  the  surface  are  large  compared  to  the  interior 
wavelength. 

Under  these  conditions,  the  wave  inside  the  coating  will  be  TEM  with 
direction  of  propagation  inward  and  normal  to  the  surface.  Hence, 

n'-£  =  n'*fi  =  0  (18) 

and  the  surface  impedance  ,  Z,  will  be  independent  of  the  polarization  of  the 
incident  radiation  and  will  be  given  by: 


Z  =  jn  tan  Npd 


(19) 


which  is  the  expression  that  would  be  obtained  for  Z  without  the  above 
assumptions  but  for  normal  incidence.  Finally,  it  should  be  noted  that  when 
dealing  with  coated  edges.  Equation  (19)  does  not  hold  at  the  edge  or  within  a 
few  skin  depths  (of  the  coating  material)  from  the  edge. 

Substituting  Equations  (18)  in  Equations  (15)  and  (16)  yields  the 
following  expressions  for  the  surface  currents: 

-ft(ft)  =  Znxt1^)  -  2nx  ^{^-(nxft)  G(r)  +  M'  G(r)}  ds ' 

r 

3($)  =  2nx?T($)  +  2nx  f  j0Zr(nx3)  G(r)  +  3xV'G(r)}  ds' 
whereas,  from  Equations  (9)  and  (10),  the  scattered  fields  are  given  by 

£S(&)  -  -/  {jen03  G(r)  +  ftxv‘  G(r)>  ds1  (22) 

Hs(j$)  =  /  {  -  —  M  G(r)  +  3xv'  G(r)J  ds’ 
s  no 

For  a  perfect  conductor,  ft  =  nx£  =  n  •  ft  =  0  and  a  similar  analysis  to  that 
given  above  leads  to 


0x^(3)  =  nx  f  { j6n  3  6(r)  -  —■  ($'  •  3)  G(r)>  ds1  (23) 

s  p 

3($)  =  2nxft’($)  +  2nx  /  3xV'  G(r)  ds'  (24) 

s 

£S($)  =  -  J(jSn  3  G(r)  -  ~  ($•  .  3)  G(r )}  ds'  (25) 

s  p 

ftS(&)  =  /  3xV'  G(r)  ds'  (26) 

s 


(20) 

(21) 
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3.  THE  AS-MM  HYBRID  DIFFRACTION  TECHNIQUE 

The  general  theory  of  the  AS-MM  hybrid  diffraction  technique  is  described 
in  detail  in  [4].  This  technique  is  based  on  the  MFIE  as  given  by  Equation 
(24).  The  basic  approach  is  to  divide  the  total  surface  S  of  the  body 
into  S*  and  Ss,  which  represent  the  lit  and  shadowed  regions,  respectively. 
The  currents  in  these  regions  are  given  by  (26),  which  can  now  be  written  as: 


3*(R)  =  2nxiH1($)  +  2nx  -f  3£x7'  G(r)  ds'  +  2nx  f  3Sx7'  G(r)  ds' 


3S($)  =  2nxfi1(i$)  +  2nx  {  3£x7'  G(r)  ds'  +  2nx  f  3Sx7‘  G(r)  ds' 


Each  region  is  further  divided  into  an  asymptotic  (AS)  region  and  a  moment 
method  (MM)  region.  The  MM  regions  are  taken  to  be  small  regions  on  either 
side  of  the  shadow  boundary  ana  the  AS  regions  as  the  remaining  regions 
in  S4  and  Ss.  In  order  to  be  able  to  solve  for  3*  and  3s  in  the  AS  and  MM 
regions,  it  is  assumed  that  3s  in  Equation  (27)  is  small  so  that  it  can  be 
neglected.  This  defines  an  approximate  lit  side  optics  current  3*p  given  by 


=  2nxfi1  (ft)  +  2  nx  j  3*xV'  G(r)  ds' 


Then  replacing  3£  in  (28)  by  3^p  defines  an  approximate  shadow  side  optics 


current 


3^  (R)  =  2nxiV(ft)  +  2nx -f  3*  xV'  G(r)  ds'  +  2nx  j  3*  x7'  G(r)  ds' 


The  exact  currents  in  Sl  and  Ss  are  then  given  by 


(31) 


34(ft)  =  3jp(i5)  +  t*(8) 

3S(^)  -  ^(S)  +  fS(*)  (32) 

Substituting  Equations  (31)  and  (32)  into  the  MFIE  (27)  and  (28)  and  using 
(29)  and  (30),  the  following  expressions  for  the  correction 

to 

and  1  are  found: 


TA(i$)  =  2nx/  3*  xV'  G(r)  ds'  +  2nx/  T^x?1  G(r)  ds'  +  2nxf  tSxV'  G(r)  ds'  (33) 
Ss  sl  Ss 


tS(ft)  =  2nx  f  lSxV'  G(r)  ds'  +  2nx  f  f*xv’  G(r)  ds1 


(34) 


The  process  that  is  defined  for  solving  the  above  system  of  equations  is  as 

follows.  is  computed  by  iterating  Equation  (29)  with  2nxfi1(R)  the  initial 

current.  Then  3Qp  is  computed  by  iterating  Equation  (30) 

using  2nxfi1(ft)  +  2nx  /  3ppxV  G(r)  ds’  as  the  initial  current.  Equations  (33) 

and  (34)  are  now  written  down  for  SMM,  S^,  S^  and  SAS-  The  correction 

currents  in  the  MM  regions  are  now  computed  by  neglecting  the  integrals  over 

the  correction  currents  in  the  AS  regions.  It  is  argued  that  without  ignoring 

o  S 

those  current  terms  containing  integrals  over  1^  and  1^  in  the  expressions 
for  1^M(ft)  and  lj^(8),  it  is  not  possible  to  derive  sets  of  equations  which 
can  be  solved  by  the  moment  method.  The  correction  currents  in  the  AS  regions 
can  now  be  computed  by  integrating  over  the  moment  method  correction  currents, 
again  ignoring  those  current  terms  containing  integrals  over  and  Ta<j. 

The  above  process  defines  a  first-order  approximation.  A  second-order 

tO  -k  $ 

and  1  in 

the  MM  and  AS  regions  in  the  integrals  for  t*M($),  t^M(ft),  T^(fi)  and  !A<-(8) 
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Higher  order 


but  this  time  without  neglecting  the  integrals  over  and 
approximations  are  obtained  in  a  similar  fashion. 

It  was  found  in  [4]  that  it  was  necessary  to  use  second-order,  and  in 
some  cases,  third-order  approximations  to  yield  a  good  value  for  the 
current.  However,  no  criterion  was  found  for  fixing  the  size  of  the  MM 


regions. 


4.  THE  ITERATION  METHOD 

The  iteration  method  is  based  on  the  AS-MM  hybrid  technique  and  is  an 
attempt  to  do  away  with  the  MM  regions  of  the  latter.  Hence,  the  exact 
current  is  written  as 

3(8)  =  3op(S)  +  1(8)  (35) 


\* 

1 

l 

r* 

K 

* 


and  the  surface  of  the  scatterer  is  divided  into  lit  and  shadow  regions 
separated  by  the  geometrical  optics  shadow  boundary.  3*p(8)  dnd  3^(8)  dre 
given  by  Equations  (29)  and  (30),  respectively,  and  are  solved  by  iteration. 

Q  c 

The  initial  currents  for  S  and  Sa  being  the  same  as  those  of  the  AS-MM  hybrid 
method. 

Having  solved  for  the  optics  currents,  the  correction  currents  are  now 
considered.  Instead  of  using  MM  regions  as  in  AS-MM,  Equations  (33)  and  (34) 
are  now  solved  by  iteration  just  as  the  optics  current  were.  Only  one 
approximation  is  called  for  here  and  that  is  the  neglecting  of  the  integral 
over  Is  in  Equation  (33).  In  second-  and  higher-order  approximations,  this 
integral  is  put  back  in  Equation  (33). 

The  full  set  of  equations  for  the  iteration  method  will  be  given  below. 
These  equations  are  based  on  Equations  (20)  and  (21)  describing  the  magnetic 
and  electric  surface  currents  induced  on  a  perfect  conductor  coated  with  a 
thin  layer  of  material  with  a  very  large  refractive  index.  The  magnetic 
correction  currents  will  be  denoted  by  8.  Hence,  the  exact  magnetic  current 
is  given  by 


W)  =  fiop(ft)  +  W) 


(36) 
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Since  the  equations  are  rather  cumbersome  for  the  general  case,  the  following 
two  operators  will  be  defined: 


L[3]  =  2nx  /  3x7'  G(r)  ds‘  (37) 

s 

K[3]  =  2nx  j-  j0(nxJ)  G(r)  ds‘  (38) 

s 

The  optics  currents  are  given  by: 

^p  =  2^  ♦  K[3'pZr]  *  L[3jp]  (39) 

5op  *  *  ^opM  *  Lt3op3  *  ^op^  *  L[3Jp]  (40) 

*'p  .  -Z^.E1  *  K[ft*op/Zr]  *  L[»‘p]  (41) 

Sop  ’  -2"«Ei  *  «fi*p/Zr]  ♦  L[«Jp]  •  K[S^p/Zr]  *  L[S‘p]  (42) 

The  correction  currents  are  given  by: 

T*  =  K[tZZr]  +  L[t4]  +  K[3^pZr]  +  L[3^p]  +  {K[tSZr]  +  L[tS]>  (43) 

tS-=  K[  TAZr]  +  L[t*]  +  K[tSZp]  +  L[tS]  (44) 

fi*  =  K[ft£/Zr]  +  L[fi‘]  +  K[ft*p/Zr]  +  L[^p]  +  ( K[f5S/Zr]  +  L[fiS]}  (45) 

=  K[^/Zp]  +  L[fi£]  +  K[?iS/Zr]  +  L[fiS] 


(46) 


The  case  of  a  perfect  conductor  can  be  retrieved  from  these  equations  by 
setting  ft  =  ft  =  0  and  then  Zr  =  0.  The  terms  in  the  braces  are  set  to  zero 
for  the  first-order  calculation. 
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5.  THE  SQUARE  CYLINDER 

The  first  application  of  the  iteration  method  to  be  considered  is  that  of 
a  two-dimensional  cylinder  of  a  square  cross  section.  The  computer  program 
for  the  iteration  method  was  set  up  to  solve  for  the  electric  currents  only, 
since,  as  can  be  seen  from  Equations  (39)  -  (46),  the  magnetic  currents  can  be 
obtained  from  the  same  program  by  simply  performing  the  duality 
transformati on: 


The  geometry  under  consideration  is  shown  in  Figure  1.  The  computer 
program  uses  the  angle  of  incidence  e1,  although  the  complementary 
angle  a  will  also  be  referred  to  here.  For  a  two-dimensional  scattering 
problem 


G(r)  ■  ir  Ho?l(Br) 


v'G(r)  •  -  !j  n|2)(sr) 


(21  (21 

where  Hv  '  and  H.v  '  are  the  Hankel  functions  of  the  second  kind  of  order 
o  1 

zero  and  one,  respectively.  Due  to  symmetry,  the  angle  of  incidence  can  be 


restricted  to 


where  the  case  of  grazing  incidence  has  been  omitted  for  the  time  being.  With 
the  above  restriction,  sides  1  and  4  are  the  lit  sides,  whereas  sides  2  and  3 
are  the  shadow  sides.  The  distance  d  along  the  perimeter  of  the  cylinder  is 
measured  in  an  anticlockwise  direction  from  the  origin  of  the  coordinate 
system  which  is  situated  at  the  corner  A. 

The  incident  field  is  given  by 


tji  =  \  e’j*u  ‘  ft 


where,  u  is  the  direction  of  propagation  of  the  incident  field  and  ft  is  the 
radius  vector  to  the  observation  point. 

For  the  case  of  a  perfectly  conducting  square  cylinder,  as  for  any 
structure  comprised  of  planar  surfaces,  the  principle  value  integrals  reduce 
to  ordinary  integrals  due  to  the  fact  that  the  L[3],  or  L[t],  terms  are 
identically  zero  for  the  integration  and  the  observation  point  on  the  same 
side  of  the  cylinder.  However,  the  introduction  of  a  finite  conductivity 
gives  rise  to  and  K[ft]  terms  for  which  the  principle  value  does  have  to 

be  taken.  The  calculation  of  these  terms  is  discussed  in  Appendix  B. 

Figures  2  and  3  show  the  magnitude  and  phase,  respectively,  of  the 

electric  current  on  a  .705A  wide  cylinder  with  a  surface  impedance 

of  Z  =  (.56,  -  .16).  The  computation  was  carried  out  for  the  first  order 
r 

only.  Due  to  the  poor  result  as  compared  to  the  method  of  moments,  MM, 
computation  (considered  exact)  [5],  the  same  run  was  made  again,  but  this  time 
for  a  perfect  conductor.  The  results  are  shown  in  Figures  4  and  5.  Despite 
the  fact  that  the  overall  agreement  is  much  better  this  time,  there  is  still 
an  unexplained  discrepancy  at  the  lit  edge  A. 


In  an  attempt  to  solve  this  specific  problem  and  at  the  same  time  to 
retain  the  problem  at  a  reasonable  level  of  complexity,  the  case  of  interest 
was  restricted  to  that  of  a  perfect  conductor. 

Clearly,  there  is  some  difficulty  involved  in  computing  the  currents 
close  to  the  edge  of  the  cylinder.  In  order  to  understand  this  problem, 
consider  what  happens  when  computing  the  optics  current  at  ft  on  side  1: 

3op(£l)  =  2nxHl  +  2nx  K  ^op^V  x^‘  G(r)  d*4  (52) 

where  ft  =  Jl^x  and  ft*  =  t^y.  Then,  for  field  points  close  to  the  edge  A 
( i .  >  0)  and  integration  points  close  to  A  («.^  +  0),  the 
distance  r  =  |ft  -  ft' |  *  0  and  the  Hankel  function  has  the  following  small 
parameter  behavior: 


h[2) ( 3r) 


2 

ti  t3  r 


(53) 


Numerical  integration  of  the  inverse  distance  function  1/r  as  r  +  »  calls  for 
very  samll  subdivisions  in  order  to  minimize  the  error.  The  closer  the  edge 
is  approached,  the  smaller  the  subintervals  would  have  to  be.  This  process  is 
clearly  unfavorable  since  the  computation  time  would  increase 
correspondingly.  An  analytical  evaluation  of  Equation  (52)  (for  details  see 
Appendix  C)  as  >  0  and  *  0  yields 

V*I  *  0)  =  2H^1  *  °)  -  7  Jop(*4  +  0)  (54) 

Numerically  integrating  the  integral  in  (52)  for  points  close  to  the  edge,  but 
avoiding  numerical  instability,  would  mean  that  the  minimum  possible 
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values  of  z,  and  z ^  would  be  about  0.02a  and  hence  the  edges  would  be 
truncated.  In  Figure  4,  the  minimum  value  of  approach  to  the  edges,  6,  was 
taken  to  be  6  =  . 0 1 X  giving  an  all  around  instability.  On  the  other  hand, 

Equation  (54)  sets  no  lower  bound  on  the  value  of  6,  so  that  the  edges  can  be 

left  intact. 

The  result  of  numerically  integrating  Equation  (52)  with  the  minimum 
values  of  and  given  by  6  =  0.02A  is  shown  in  Figures  6  and  7  which  show, 
respectively,  the  magnitude  and  the  phase  of  the  current  induced  on 
a  0.6A  side  square  cylinder  for  a  =  45°  .  Notice  the  overall  discrepancy 
between  the  current  computed  by  MM  and  that  computed  by  the  iteration 

method.  Figures  8  and  9  show  the  vast  improvement  obtained  by  using  Equation 
(54)  for  the  iteration  method,  *'ii,  the  minimum  values  of  and  z ^  taken  to 
be  6  =  10~^A.  This  value  was  henceforth  taken  as  the  fixed  minimum  value 
for  6.  The  influence  of  the  order  of  the  approximation  can  be  seen  by 

comparing  Figures  10  and  11,  a  first  order  calculation,  with  Figures  8  and  9, 

a  third  order  calculation.  As  can  be  seen,  there  is  an  overall  improvement  in 

the  value  of  the  magnitude  of  the  current  in  going  from  first  order  to  third 

order,  except  in  the  deep  shadow  region. 

Further  third  order  results  for  5  =  10~6A  are  shown  in  Figures  12  to 

21.  Figures  12  to  15  are  for  a  =  46°,  w  =  2.4a  and  w  =  3.7a,  whereas  Figures 
16  to  21  are  for  a  =  65° ,  w  =  .6X,  2.4A  and  3.7a. 
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6. 


IMPROVING  THE  ITERATION  METHOD 


The  final  result  of  any  iteration  method  is  dependent  on  the  initial 

value  used  to  start  off  the  iteration  procedure.  As  already  explained  in 

Chapters  2  and  3,  the  present  iteration  method  is  based  on  the  AS-MM  so 

that  is  iterated  first  with  an  initial  value  of  2nx?P  and  then  3s  is 
op  op 

iterated  with  an  initial  value  of  2nxiV  +  2nx  /  3gpX^  G(r)  ds'.  Figures  22 
and  23  show  the  resulting  magnitude  of  the  optics  current  as  compared  with  the 
MM  for  W  =  ,6A  and  3 . 7 A  and  a  =  65°.  Note  that  in  both  cases,  the  shadow  side 

optics  current  is  much  closer  to  the  true  current  than  the  lit  side  optics 

current. 

The  process  described  above  involves  taking  the  shadow  side  current  as 

zero  when  iterating  the  lit  side  current.  This  explains  the  relatively  smooth 

slow  varying  result  obtained  for  the  lit  side  optics  current,  for  there  are  no 

iterations  with  the  shadow  side  at  this  stage.  On  the  other  hand,  when 

computing  the  shadow  side  current,  interaction  with  the  lit  side  is  taken  into 

account.  This  explains  why  the  shadow  side  optics  current  follows  more 

closely  the  true  current  than  does  the  lit  side  optics  current. 

A  possible  solution  to  this  problem  would  be  to  start  off  by 

computing  3*  and  taking  3*  as  20XH1.  This  would  allow  for  interaction 

between  the  shadow  side  and  a  reasonable  approximation  for  the  lit  side. 

Having  iterated  3*  ,  the  lit  side  optics  current  can  be  computed 

using  3*  which  would  give  a  much  better  result  than  formerly  when  3s  was 

op  J  op 

taken  as  zero.  This  approach,  as  well  as  others,  will  be  tried  out  in  the 


7.  CONCLUSION 


An  iteration  method  for  solving  the  MF I E  and  the  EFIE  for  a  general 
scatterer  has  been  presented.  This  method  is  similar  to  the  AS-MM  hybrid 

technique  in  that  the  true  surface  current  is  written  as  the  sum  of 

approximate  optics  currents  and  correction  currents  on  the  lit  and  shadow  side 
regions  of  the  body.  The  two  methods  differ  in  the  way  in  which  the 

correction  currents  are  computed.  In  the  hybrid  technique  the  lit  and  shadow 

areas  are  further  subdivided  into  moment  method  and  asymptotic  regions,  the 

former  being  the  smaller.  The  MM  region  correction  currents  are  then  solved 
by  using  the  method  of  moments,  while  neglecting  their  interaction  with  the 
asymptotic  regions.  The  asymptotic  correction  currents  are  then  obtained  by 
integrating  the  moment  method  region  currents  numerically,  while  again 

neglecting  the  interaction  between  the  asymptotic  region.  These  interactions 
are  accounted  for  in  second-order  and  higher  approximations.  In  the  iteration 
method,  the  correction  currents  are  obtained  by  iteration  just  as  ♦  optics 
currents  are.  The  only  approximation  made  when  computing  the  correction 
currents  via  the  iteration  method  being  the  neglecting  of  the  contribution  of 
the  shadow  side  correction  current  when  computing  the  lit  side  correction 
current.  Again,  in  second-  and  higher-order  approximations  this  contribution 
is  taken  into  account. 

As  an  example,  the  case  of  a  two-dimensional  cylinder  of  square  cross 
section  was  considered.  The  currents  obtained  by  the  iteration  method  were 

compared  to  the  currents  obtained  by  the  method  of  moments.  The  overall 

excellent  agreement,  except  for  small  regions  in  the  shadow  region  for  large 
bodies,  suggests  that  the  order  of  the  iteration  process  should  be  changed, 
viz.,  the  shadow  side  optics  currents  should  be  computed  first  and  then  the 
lit  side  optics  currents.  This  small  discrepancy  in  the  shadow  region  does 


not,  however,  minimize  the  importance  of  the  iteration  method  which  has 
successfully  done  away  with  the  necessity  of  a  moment  method  region  on  the 
surface  of  the  scatterer. 

The  formalism  was  shown  to  be  adaptable  to  problems  of  non-perfect 
conductors,  and  the  special  case  of  a  coated  perfect  conductor  was  considered. 

Finally,  it  is  pointed  out  that  the  iteration  method  does  not  give 
erroneous  results  for  geometries  having  an  interior  resonance  due  to  the  fact 
that  the  inhomogeneous  integral  equation  for  the  surface  current  is  formulated 
as  two  inhomogeneous  integral  equations,  one  for  the  optics  currents  and  one 
for  the  correction  currents  and  neither  of  these  currents  alone  can  sustain  a 
zero  tangential  field  at  the  surface  of  the  scatterer.  The  iteration  method 
has  another  advantage  over  the  moment  method  in  that  it  can  be  used  for 
electrically  large  bodies  that  the  MM  cannot  handle,  since  it  is  not  limited 


by  computer  storage. 
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APPENDIX  A 


DEFINITION  OF  ELECTROMAGNETIC  QUANTITIES 

The  complex  dielectric  constant  is  given  by 


je"  =  e  (e‘  -  je")  =  e  e 
u  o'r  *  r'  or 


j  tan6  ) 

e 


The  complex  magnetic  permeability  is  given  by 


u  =  u*  -  ju"  =  uq ( u -  jpp  =  MQur  =  J.(  1  -  j  tan6^) 

A  prime  denotes  the  real  part  of  a  complex  quantity 
A  double  prime  denotes  the  imaginary  part  of  a  complex  quantity 
A  small  subscript  o  refers  to  vacuum 

A  small  subscript  r  means  that  quantity  relative  to  the  same  quantity 
vacuum. 

e 1  =  — -  =  dielectric  constant 
r  £o 

e"  =  —2—,  where  o  =  conductivity 
r  e  a)  J 

0 

Up  =  —  =  relative  permeability 
uo 

n 

Up  =  ■— —  =  relative  magnetic  loss  factor 
6^  is  the  magnetic  loss  angle 
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APENDIX  B 


CALCULATION  OF  THE  PRINCIPLE  VALUE  INTEGRALS 
The  integrals  defined  by  Equation  (38)  are  singular  as  r  -#•  o.  These 
integrals  are  computed  by  dividing  them  into  subintervals  which  are  integrated 
separately.  The  subinterval  containing  the  singularity  is  of  the  form 

6/2  /p\ 

I  =  J  (M)  d£  (B.l) 

-6/2  0 

and  is  computed  by  taking  its  principle  value.  Here,  £  is  a  length  parameter 
and  5  is  the  size  of  the  subinterval. 

The  accuracy  of  the  above  integration  depends  on  the  width,  6,  and  on  the 
approximation  used  for  the  asymptotic  expansion  of  the  Hankel  function.  The 
usual  small  argument  expression  used  for  H^  is  [6]: 

H^2)(0£)  ^1  -  j|  zn(|  y0£)  (8.2) 

substituting  (B.2)  in  (8.1)  yields  a  principle  value  of 

1^6(1  -  j~  fcn(^~)}  (B.3) 

where  £ny  =  .57721566  is  Euler's  constant  and  e  is  the  natural  base  of 
logarithms. 

The  values  of  as  given  by  (B.2)  can  be  compared  with  its  exact 

values  and  it  is  found  that  Eq.  (8.2)  is  in  error  for  (£/A)  >  .02,  giving  rise 
to  an  error  in  Eq.  (8.3).  From  asymptotic  analysis  [7],  we  know  that  we  can 
use  larger  values  of  (£/A)  and  decrease  the  error  by  including  higher  order 
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terms.  The  asymptotic  expansion  for  '  including  terms  of  the 
order  (jit)2  is  [8] 

H^2)(3)  «te[l  -  -j(fU)2]  Cl  -  j|  tn(|  Y0t)]  -  ^r(3*)2  (B.4) 

yielding  the  following  expression  for  the  principle  value  of  the  integral  I: 


I  -  io(e«)2]  Cl  -  K  *"(4=^)3  -  7?t(66)2} 


(B.5) 


Equation  (B.4)  is  accurate  for  values  of  (t/A)  <  .1.  For  a  value  of  6  =  .12A, 
Eq.  (B.3)  has  an  error  of  1.4%  in  comparison  to  Eq.  (5),  which  can  be 
considered  exact  for  this  value  of  6. 
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APPEHDIX  C 


CALCULATION  OF  THE  EDGE  CURRENTS 


In  this  appendix,  the  result  given  in  Eq.  (63)  will  be  derived.  The 
derivation  will  be  given  for  a  wedge  of  internal  angle  a.  The  geometry 
involved  and  the  definitions  of  the  various  quantities  are  described  in  Figure 
24.  The  equation  to  be  considered  is 


3A(£)  =  2nx?i’(£)  +  2nx  /  3B(£‘)x$'  G(r) 


(C.l) 


Now  consider  the  integral 


t(£)  =  2nx  /  <5 B  ( £ ' )  xv'  G(r)  d£' 


(C.2) 


where 


3B(£')  =  JB(£‘)t',  G  =  H^2) (tsr)  r' 


therefore 


1(£)  =  y~  /  nx(t'xr')JB(£')H^(gr)  d£' 


(C.3) 


nx(t'xr')  =  nx[t'x{(r‘  *  t')t'  +  (r'  •  n')n'}] 


=  t(r'  •  n' )  =  t  cosy 
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T(t)  =  4^  t  /  JB(£ • )  H.^COS*  d£' 

^  n 1 !■ 


(C.4) 


where  cos^  =  —  sina 


(C.5) 


/  2  2 

and  r  =  /  V  +  £  -  2«.  £'cosa 


(C.6) 


Using  the  small  argument  expansion  for  H' 


we  get. 


t<o  -  &i  i  At 


' )  j  —1 
1  J  ngr 


£  Sincx 
0 


where  L  is  much  greater  than  t  and  gr  is  much  smaller  than  unity.  The  above 
expansion  for  the  Hankel  function  can  be  used  since  we  are  only  interested  in 
a  small  neighborhood  around  the  edge.  Therefore, 


t(D  - 


-  t  J8(£*  *  0)  ^  j  - - 

11  0  +  %C  -  2£'£C0Scx 


(C.7) 


This  is  a  standard  integral  which  can  be  integrated  to  give 


7  A*'  ♦  o)  tan'1(~£sf^a)  -  tan_1(-  cote.) 
- 


therefore. 
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Figure  3.  Phase  of  current  on  cylinder 

W  =  0.705A ,  a  =  45°  and  Z  =  (.56,  -.16) 

_  MM 

---------  Iteration  Method 
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Phase  of  current  on  cylinder 
with  edge  truncation,  6  =  0.02A 
W  =  0.6A ,  a  =  45° 
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Figure  8.  Magnitude  of  current  on  cylinder 

W  =  0.6A,  a  =  45° 

_  MM 

--------  Iteration  Method  (3rd  order) 
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Figure  16. 


Magnitude  of  current  on  cylinder 
W  =  0.6X ,  a  =  65° 

_  MM 

---------  Iteration  Method 
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Figure  17. 


Phase  of  current  on  cylinder 
W  =  0.6A,  a  =  65° 

MM 

.  ..  Iteration  Method 


Figure  20. 


Magnitude  of  current  on  cylinder 
W  =  3.7\,  a  =  65° 
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_________  iteration  Method 
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